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Abstract 


We evaluate in detail the production of deuterium, tritium, 
and helium-3 from nuclear reactions of accelerated charged particles 
with the ambient solar atmosphere. We use updated cross sections 
and kinematics, we extend our calculations to very low energies 
(~ 0.1 MeV/nucleon) , and we calculate the angular distribution 
of the secondary particles. We then compare the calculations with 
data on accelerated isotopes from solar flares. In particular, we 
consider the August 1972 events for which both He'* and nuclear gamma- 
rays were observed. Finally, we provide an explanation for He 3 — rich 
events in terms of the angular distributions of secondary isotopes, 
and we make predictions on the flux of 2.2 MeV gamma rays from such 


flares. 
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I . Introduction 

The isotopes H 2 , H 3 and He 3 were detected in solar cosmic rays 
with abundances much larger than the H 2 and He 3 abundances or upper 
limits in both the solar wind and the photosphere (for summary of 
observational data see Section IV). Therefore, it is generally 
believed that energetic H 2 , H 3 and He 3 from solar flares are produced 
mainly by nuclear reactions of accelerated H 1 and He 4 with the ambient 
material in the flare region. The principal difficulty with this 
interpretation, however, appears to be caused by observations of 
flare associated events with very large He 3 /He 4 ratios, but which do 
not show detectable H 2 and H 3 fluxes (Garrard, Stone and Vogt 1973, 

Anglin 1974) . In order to understand this apparent discrepancy in the 
context of all the available information on nuclear reactions in flares 
(both isotopic and gamma-ray data) , we have reexamined in detail the 
production of H 2 , H 3 and He 3 by accelerated particles, and we have 
compared our results with these data. 

Previous detailed calculations on H and He isotope production and 
applications to galactic cosmic rays were done by Ramaty and Lingenfelter 
(1969) and by Meyer (1971, 1972). In the present paper we apply these 
calculations to solar flares. We also extend their treatments of the 
nuclear physics as follows: We calculate the secondary spectra down to 

very low energies (~ 0.1 MeV/nucleon) , and we derive angular distribu- 
tions for the isotopes; we supplement the cross sectional data on pCNO 
reactions at low energies; and we consider the destruction of the 
isotopes (in particular tritium) by nuclear reactions. The contributions 
of pCNO reactions, however, turns out to be unimportant for most cases 


of interest. 
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Isotope production in solar flares without consideration of reaction 
kinematics was previously treated by Lingenfelter and Ramaty (1967) . 

We then compare our results with data. For the flares of August 
1972 for which both He 3 and nuclear gamma rays were observed (Webber 
1974; Chupp et al. 1973), the isotope data implies that the product of 
the interaction time of the particles and the density of the ambient 
medium was about 4xl0 13 cm -3 £; this corresponds to a matter traversal 
for relativistic particles of ~ 2 gem 3 . A similar matter traversal 
(~ 1 gem -3 ) is required to account for the deuterium and tritium obser- 
vations from other flares. The He 3 observations from He 3 -rich flares, 
however, require a much larger path length ( >1 0 gem 3 ) ; moreover, because 
of the absence of H 3 and H 3 , these events cannot be understood in terms 
of isotropic nuclear interactions in flares. We show that if the primary 
particles are directed downwards toward the Sun, the depletion of H 3 and 
H 3 could result from the preferential emission of these isotopes into 
the forward hemisphere, and from the destruction of H 3 by nuclear 
interactions. We expect such anisotropies mainly for small flares. 

The validity of our model for He 3 -rich flares could be checked by 
gamma-ray observations. Our prediction is that these flares should 
produce a time-integrated 2.2 MeV flux of about 3 photons cm -3 . We 
note, however, this value is quite uncertain, mainly because of uncer- 
tainties in the evaluation of the total number of He 3 nuclei released 


from the flare. 
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II. Isotope Production 

The principal deuteron, triton and helium-3 producing reactions in 
proton-hydrogen (pp) , proton-helium (pa) and alpha-particle-hydrogen 
(o?p) Interactions are listed In Tahle 1. We also give the reaction 
thresholds and minimum secondary particle energies, E g (min) (Ramaty 
and Lingenfelter 1969) . In addition to these reactions we consider 
interactions of accelerated protons with ambient CNO, and the inter- 
actions of accelerated CNO particles with ambient hydrogen. We have 
neglected ata reactions, because no reliable data on cross sections 
exists for these reactions. We note, however, that by using approxi- 
mate cross sections, Ramaty and Lingenfelter (1969) found that aa 
reactions do not contribute appreciably to the total H 3 and He 3 pro- 
duction by galactic cosmic rays. In solar flares the contribution 
of aa reactions will be even less, because here the ratio of accelerated 

alpha particles to accelerated protons is smaller than in the cosmic 
r ays . 

For the composition of the ambient solar atmosphere we use 
H:He:C:N:0 = 1:0.07: 3.7xl0 -4 : 1.2xl0“ 4 : 6.8xl0“ 4 (Cameron 1973). For 
the accelerated particle populations we assume both power-law and 
exponential spectra, i.e. 

V E > - k i E_s (1) 

and 

N^P) = k i exp(-P/Po) , ( 2 ) 

respectively. Here N^(E) and N^(P) are the number of accelerated 
particles of kind i in the flare region per unit energy per nucleon, 

F, or unit rigidity, P; the k^s are constants determined by normalizing 
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the N^'sto 1 proton of energy greater than 30 MeV, and by using 
the composition of the ambient solar atmosphere; and s and Po are, 
respectively, the spectral index and characteristic rigidity assumed 
to be the same for all accelerated particle components. 

For isotope production in pa reactions we use the cross sections 
of Meyer (1972). For H 2 and H 3 production in pCNO reactions, and 
for He 3 production in pCNO reactions above about 20 MeV/nucleon we use 
the cross sections of Ramaty and Lingenfelter (1969) . The production of 
He 3 below 20 MeV/nucleon comes from the reaction N 14 (p, He 3 )C 12 which 
has a very low threshold energy (4.71 MeV) . At 13 MeV the data of 
MacLeod and Reid (1966) gives a cross section of 28 mb. At lower 
energies we calculate the cross section from data on the inverse reaction, 
C 12 + He 3 -> N 14 + p (Kuan, Bonner and Risser 1964) by using the detailed 
balance principle. After averaging over resonances, some representative 
values are: a(7MeV) =*■ 1.9mb, a(8MeV) ^ limb, a(lOMeV) ^ 18mb . Above 
25 MeV/nucleon the data merges with the cross sections given by Ramaty 
and Lingenfelter (1969) , 

Beside the production cross sections, we also need the angular and 
energy distributions of the secondary products. For the 2-body reaction, 
p+a hH 2 + He 3 , we use a fit to the angular distributions as given by 
Ramaty and Lingenfelter (1969). An important feature of this distri- 
bution is that in the c.m. (center-of-mass) frame the deuteron is emitted 
in the forward direction, while the He 3 is coming out in the backward 
direction. This anisotropy becomes more pronounced with increasing 


projectile energy. 
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For the rest of the pc* reactions, where there are 3 or more parti- 
cles in the final state (3-body reactions )> the available experimental 
data is less complete- Ramaty and Lingenfelter (1969) assumed that 
for these reactions the velocity of the secondary particle in the c-m- 
frame is always less than or equal to the velocity of the center of 
mass. As a consequence * in the laboratory frame the secondary particle 
is always emitted in the forward direction- This can be compared with 
the data: At an incident energy of 300 MeV the angular distribution 

of deutrons from 3-body reactions is strongly peaked in the forward 
direction in the laboratory frame (Innes 1957)- At 95 MeV, the He 3 
angular distribution is also peaked in the forward direction (Tannenwald 
1958) , but at 300 MeV, a substantial fraction of the He 3 nuclei are 
emitted in the backward direction as well (Innes 1957). Thus, the 
assumption of Ramaty and Lingenfelter (1969) gives correct results 
regarding the concentration of deuterons in the forward direction, but 
it underestimates the amount of He 3 produced in the backward direction. 
There is no data on the angular distribution of tritium, but in analogy 
to He 3 , we expect that the above assumption would also underestimate 
the amount of tritium emitted in the backward direction. 

We perform our calculations for two limiting models: a thick- 

target model in which the accelerated particles move from the flare 
region downwards into the Sun, undergoing nuclear interactions as they 
slow down in the solar atmosphere; and a thin- target model in which the 
spectrum of accelerated particles is not modified during the time in 
which the nuclear interactions take place- The latter model assumes 
that either the total path length traversed by the particles at the 
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Sun is small in comparison with their interaction length, or that the 
particle energy loss from ionization and nuclear interactions is just 
balanced by energy gains from acceleration. 

Let us first consider the thick-target model. For this model we 
assume that the primary particles are collimated into a beam directed 
downward into the Sun. We define the direction of the beam as the 
forward direction. Because we can assume that the primary particles 
are slowed down mainly by ionization losses, the total number of 
secondary particles per unit energy per nucleon and unit solid angle 
produced from a particular type of reaction during the slowing down 
of the primary particle is given by 


dN s (E g> p) _ nt f dE " N <E / )r E i dx d s a(E,E*,p*) f 3(E s ,p) l ' 1 (3) 

dQ n H J 0 4 dE (E) dE* dtf U(E?,p*)J dE ‘ 

s ® 

A A 

Here (E g ,p) and (E g ,p ) are the energies per nucleon and directional 
cosines of the secondary particle in the frame of the Sun and the c.m. 
frame, respectively, with p = 1 and p = 1 corresponding to the forward 
direction; n t and n^ are the target nuclei and proton densities in the 
ambient medium and mp is the proton mass; N(E) is the initial number of 
primary particles in the beam per unit energy per nucleon; dE/dx is the 


stopping power for the primary particles in the solar atmosphere; 

ml ^ 

d 2 a(E,E , p ) /dE dQ Is the differential cross section in the c.m. frame 
s s 

for producing a secondary particle of energy E g and directional cosine 


p ; and 3 (E g , p) /a (E g , p ) is the Jacobian of transformation between 
(£ s ,p) and (E*,p*). This transformation is 





* * 
P s ^ 


( 4 ) 
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n - ((yJ-D 


1/2 


* * *. . 

Y s + Vs^ )/p s» 


(5) 


where y^ is the Lorentz factor of the center of mass in the frame of 
the Sun, y s - l+E s /m p , y g = l+E s /m p , and p g and p g are momenta 
corresponding to E s and E , respectively. The Lorentz factor y c > is a 
unique function. of the incident energy E, and is given by 

V C ■ (E cm S + V - W 

where E cm is the total available energy in the c.m. frame, and 

and m t are the masses of the projectile and target, respectively. 

In evaluating equation (3) we distinguish between cases where 

the final state consists of 2 particles and cases where the final state 

& 

consists of 3 or more particles. For a 2-particle final state, E g = 

ft ft 

E m ; E m is a function of the incident energy E, 


V„ - < E L + "! - «r S >/<*»» E cm> • 


(7) 


ft ft 

where Y m = 1 + E m /m g , m g is the mass of the secondary particle under 
consideration, and m r is the mass of the other secondary particle. For 
a 3-particle final state, E g has a range of values up a maximum given 
by equation (7), where now m r is the sum of the masses of all the other 
secondary particles. 

In the absence of detailed experimental data on the angular and 

energy distributions in 3-body reactions, following Ramaty and 

Lingenfeleter (1969), we make the assumption that E* is a unique function 

s 

of E. For pa reactions this function is given by, E* = E* for, y*<Y c and 


” m s (Y c -D otherwise. For op reactions, for a given E ,E* is the 


s m 

ft 

'cm* “s 


same as the E g for pa reactions with the same E cm . 

If the c.m. energy of the secondary particle is a unique function 
of the incident energy, E, then the center-of -mass differential cross 
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section be written as 

d 3 c(E,E s ,M,*) b 1 a(E) p( E> ^*) 6(E*-E*). 


( 8 ) 


dEg dff 


2n 


Here cr(E) is the total cross section as a function of energy, and 
P(E, (a,*) dp* is the probability that in the c.m. frame the secondary 
particle will have directional cosine in dp, around p . By substi- 
tuting equation (8) into equation (3), and by interchanging the orders 
of integration, we obtain 
dN c 


k (E s , p)= ^t ct(e) dx dE PCE 

2tt 


,P*) f 9 (Eg, p) 1 _1 

n U(E*,p*)_ 


n(>e) 


(9) 


dfi n H m p dE d E^ 

where N(>E) is the initial number of primary particles in the beam 

with energies per nucleon greater than E. All quantities on the r.h.s. 
of equation (9) are functions of E and p*. For a given E g and p, E 
is obtained by solving the equation 


Y s = Y c Y £ 


-(y/ -l ) 1 / 2 


pgM'j 


( 10 ) 


where both v* and v are functions of E; the directional cosine fj, is 
s c 


thei obtained by evaluating the equation 

P* - (Y 0 P s l* -Y s (Y=-l) l/= ) /Pi 
The Jacobian in equation (9) is given by 

-* .* 


3(Eg,p) 


0 + 


d Y c P s 


*1 

J Ps 


(ID 


( 12 ) 


a<E*,p*) L W s jr 

We have evaluated equation (9) numerically for the various reactions 

listed in Table 1, and with a and P as given by the experimental data 

discussed above. These calculations were carried out for a variety of 

values for E and p, as well as for various primary spectra. Some 
s 

of these results will be presented in the discussion below. 
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If we Integrate equation (9) over dQ we obtain the total number 
of secondary patticles per unit energy interval 


N S (E S ) - 2n 


'dp, dN s(|s.^) 


( 13 ) 


For a fixed E g , the transformations of equations (4) and (5) imply that 
p is a unique function of E*, which in turn is a unique function of E 
by our assumptions. We can thus transform the integration variable in 
equation (13) from p to E. For the derivative dp/dE we find 


d ,P = 9 (E s , (1) dE* 
dE 15 TeJ77I*) 



,-l 


(14) 


Then, by substituting equations (9) and (14) into equation (13) we obtain 

W - <V”h) J" 4®) ^ dE, (15) 

0 * 

where 


f(E,E s ) = P(E,p*)(|£|)= P(E, p*) (y^-l)" i/2 (Y* m -i (16) 

is the probability that in the laboratory frame the energy of the 
secondary particle will be in dE g around Eg. Equation (15) could have 
been obtained directly from elementary consideration, but we need the 
f°™alism in order to calculate the angular distributions. 

For the thin-target model, the equations analogous to equations 
(3) and (15) are 

/p P _ r\ / n % ^ 


and 


r N(E) ± d a a(E,E s ,/) f 8(6^) 1 

d " nH H "p iTnfT-] 


dE 


(17) 


W = 


~ x x J N(E) p(E) f (E,E s )dE, 


( 18 ) 


where the total path length traversed by relativistic 
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primary particles in the interaction time t L , 3 is velocity in units 
of c, and N(E) is the number of primary particles per unit energy per 
nucleon in the interaction region. For this model, however, we do not 
need the equivalent to equation (9) , because we shall evaluate isotope 
production in the backward hemisphere for the thick-target model only. 
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III. Numerical Results 

In Figures 1, 2, and 3 we show the H s , H 3 and He 3 production 
as evaluated from equation (9) and (13) for the thick-target model. 

Here the primary particles have exponential rigidity spectra with P Q = 

150 MV. The numerals on the various curves correspond to the reactions 
listed in Table 1 and the total curve is the sum of these partial 
production modes. Except for the backward spectra, all curves in 
Figures 1, 2 and 3 represent isotope productions integrated over all 
directional cosines, -1 s p s 1. The backward spectra are obtained 
by integrating equation (13) over -1 £ p ^ 0 only. 

We first consider the spectra integrated over all angles. We see 

that isotope production below about 10 MeV/nucleon is due entirely to 

po reactions. This is a model independent result, as can be seen from 

Table 1. Here, the minimum energies of secondary particles, E (min), 

s 

for op reactions are about 10 MeV/nucleon and therefore these reactions 
cannot contribute to isotope production at lower energies. 

Above about 20 MeV/nucleon op reactions make a major contribution 
to the total H 3 and He 3 production. For deuterium the po reactions 
remain important in this energy region as well. Moreover, around 100 
MeV/nucleon most of the H 3 comes from the reaction p4p-*d+ri + It should 
be noted, however, that the relative contribution of the reactions 
depends on the assumed ratio of the flare-accelerated alpha particles to 
protons; it is known that this ratio can be energy dependent and time 
variable. It is important to realize however, that the calculated spectra 
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below about 10 MeV/nucleon are Independent of this variability. 

We now consider the backward spectra. Because of our assump- 
tions on the c.m. energy distributions, we obtain secondary particles 
in the backward direction only from reactions (1) and (10) . For 
deuterium, this result is consistent with the experimental data 
discussed above, but for H 3 and He 3 we might expect a contribution 
to the backward spectra also from the 3-body reactions. The forward 
peaking of the H 3 in the 2-body reactions is also responsible for the 
relatively low contribution of reaction (1) to the total H 3 produc- 
tion at low energies and its relatively high contribution at larger 
energies . 

In Figures 4 and 5 we present T(He 3 /H 3 ) , the ratio of the backward He 5 
to H 2 at the same energy per nucleon. These results are for the thick- 
target model, with exponential and power-law spectra, respectively. 

From the above discussion, these ratios are lower limits, because of 
the additional backward contributions to He 3 from 3-body reactions. If 
we assume that Innes 1 (1957) data at 300 MeV is representative, then about 
1/3 of the total He 3 in 3-body reactions is emitted in the backward 
direction. Then from Figure 3 we estimate that the backward He 3 /H 2 
ratios in Figures 4 and 5 could be larger by factors of 2 to 3. Another 
source of uncertainty in these ratios is caused by poor cross sectional 
data for the reaction p+a-*H 2 +lfe 3 above ~ 400 MeV. Because we have 
assumed that in this energy region the cross section is zero (Meyer 
1972), any finite value for it will further increase the He 3 /H 2 ratios 
in Figures 4 and 5. This increase, however, should be small, because 
we do not expect many flare accelerated protons at high energies. 
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In Figures 6, 7 and 8 we show the total H 2 , H 3 and He 3 production 
for the thin-target model for various spectra of the primary particles. 
These figures have been obtained from equation (18), with N(E) given 
either by equation (2) with P Q - 50MV and 300 MV, or by equation (1) 
with s = 2 and 4, and for x^ = 1 gem -3 . As can be seen, below about 
10 MeV/nucleon the secondary spectra are only weakly dependent on the 
spectra of the primary particles. This results from the fact that 
these secondary nuclei are produced mainly by 3 ‘body pa reactions . 

From the experimental data (e.g. Ramaty and Lingenfelter 1969) it 
follows that in these reactions the energy spectra of the secondary 
particles are independent of the projectile energy. At higher energies, 
however, the secondary specta tend to resemble the spectra of the primary 
particles . 
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IV. Comparison with Observational Data 

We start the comparison of our calculations with observational 
data for the August 1972 events for which both accelerated He 3 and 
nuclear gamma rays were observed. According to Webber et al. (1974) 
r(He 3 /He 4 ) , the He a /He 4 ratio at the same energy per nucleon, is about 
0.02 at 60 MeV/nucleon, and it decreases monotonically with decreasing 
energy to about 0.002 at ~ 7 MeV/nucleon. 

Using our calculations as presented in Figure 8 (together with 
additional unplotted calculations) , we plot in Figure 9 F(He /He ) 
for the thin-target model with power law primary spectra characterized 
by various values of s. At ~ 60 MeV/nucleon, r(He 3 /He 4 ) “ 0.01 xj • 

By comparing this with the observed ratio, we get that x l <=- 2 gem -3 , 
and x(60 MeV/nucleon) = p x, ^0.7 gem" 3 . Because the stopping range 
of a 60 MeV/nucleon He 3 in hydrogen is about 1 gem -2 (Barkas and Berger 
1964) , the He 3 nuclei will escape from the production region without 
much attenuation. 

Let us compare this result with the gamma- ray observations. 

According to Chupp et al. (1973), the average flux in the 2.2 MeV 
line, over a time interval which coincides approximately with the 
flash phase of August 4, 1972 flare (0623 to 0633 U.T.), was about 
0.3 photons cm -3 s -1 . The 2.2 MeV line is due to neutron capture on 
hydrogen, the neutrons being produced by nuclear interactions of 
accelerated changed particles (e.g. Ramaty and Lingenfelter 1973). We 
can estimate the average number of protons in the flare region as 
follows: From the calculations of Wang and Ramaty (1974), the efficiency 

of gamma-ray production by neutron capture on hydrogen in the solar 
atmosphere is about 0.1. The observed 2.2 MeV flux then implies an 
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average rate of neutron production of ~ 10 28 neutrons s^ 1 . Using 
available calculations of neutron production (e.g. Ramaty, Cline 
and Fisk 1972), this rate then gives that n H N^(>30MeV) 1.5 x 10 44 cm~ 3 , 
where Np(>30MeV) is the average number of protons present in the flare 
region during the gamma-ray observations. We note that, whereas the 
isotope data determines the product of the ambient density and the 
interaction time, (n^t^ =* 4 x 10 13 cm~ 3 s), the gamraa-ray data gives 
the product of the same ambient density and the average number of 
protons in the interaction region. If we combine these 2 products, we 
get for the August 4, 1972 flare 

Np (>30MeV ) /ti =“ 4 xlO 30 protons (19) 

The quantity given by equation (19) in the rate of proton release from 
the flare region. The total number of protons released is the product 
of this rate and the acceleration time, T. For the August 4, 1972 flare T 
is — 1 0 3 seconds as indicated by the microwave data (Toyokawa Observatory, 
private communication, 1972) . The total number of protons released 
above 30MeV is therefore about 4 x 10 33 . 

This number should be compared with the total number of protons 
in the interplanetary medium as deduced from charged particle obser- 
vations. From the estimates of Ramaty and Lingenfelter (1973) this quan- 
tity is about 2 x 10 34 protons with energies greater than 30 MeV. The 
discrepancy of about an order of magnitude between this result and that 
obtained above from the nuclear reactions is probably due to uncertain- 
ties in the measurement of the extremely large particle flux in the 
interplanetary medium and the estimate of the volume that contains these 
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particles. We note that the number 2 x 10 34 is based on the measure- 
ments of Kohl, Bostrom and Williams (1973). Because the proton flux 
given by these authors could have been too large by at least a factor 
of 3 (M. Van Hollebeke, private communication 1972), the above discrepancy 
could be reduced by at least this factor. 

In addition to the August 1972 events, H and He isotopes were 
detected also from other solar flares (Schaeffer and Zahringer 1962; 

Hsieh and Simpson 1970; Dietrich 1973; Anglin, Dietrich and Simpson 
1973a; Garrard, Stone and Vogt 1973). There are no reported observations 
of H 2 and H 3 from individual flares; Anglin, Dietrich and Simpson 
(1973a) only present H 2 /H 1 and H^H 1 ratios, averaged over several solar 
flares. At ~ lOMeV/nucleon, F(H 2 /H 1 ) ^ 8 x 10 5 and r(H 3 /H 1 ) =* 2 x 10 5 . 
Using our calculations as presented in Figures 6 and 7 (together with 
additional unplotted calculations) , we give in Figures 10 and 11 
r(H 2 /H 1 ) and TO^/H 1 ) for the thin-target model with power law primary 
spectra characterized by various values of For 2 <; s £ 3, the obser- 
vation r(H 3 /H 1 ) =- 8 x 10 “ 5 at 10 MeV/nucleon implies that Xj is between 
.5 to 1.5 gem” 2 , and that x(10 MeV/nucleon) - px 1 is between 0.08 to 0.2 
gem” 2 . These path lengths are consistent with the H 3 observations, 
as can be seen by using the results of Figure 11. They are also quite 
comparable with the path length deduced from the He 3 /He 4 ratio of the 
August 1972 events. We note that because the stopping ranges of 10 MeV/ 
nucleon H 3 and H 3 in hydrogen are about 0.1 gem” 2 and 0.15 gcm~ 2 , these 
nuclei can escape from the production region without much attenuation. 
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However, the average He 3 /He 4 ratio of about 2 x 10 ” 2 at ~ 10 MeV/ 
nucleon used by Anglin, Dietrich and Simpson (1973a) to compare 
with their H 2 /H 1 ratios is inconsistent with these values of . 

Indeed as can be seen from Figure 9, T(He 3 /He 4 ) =* 2 x 10 “ 2 requires 
an xj ranging from about 20 gcm“ 3 to 60 gem” 2 for 2 £ s £ 3. This 
could indicate that the average He 3 /He 4 ratios are greatly modified 
by a subset of flares which are enriched in He 3 * 

There are observations of He 3 -rich events. Thus, T(He 3 /He 4 ) 
for the flare of October 14, 1969 was about 0.3 from a few MeV/nucleon to 
about 20 MeV/nucleon (Garrard, Stone and Vogt 1973, Anglin 1974). 
Furthermore, for the flare of July 30, 1970 the He 3 /He 4 ratio was about 
0.5 from about 10 to 20 MeV/nucleon (Anglin 1974). A common feature 
of these He 3 rich events is absence of H 3 and H 3 . The lower limits 
on f*(He 3 /H 3 ) and r(He 3 /H 3 ) are, in some cases, as large as 20. 

We wish to suggest that in these He 3 -rich events we observe pre- 
dominantly the products of nuclear interactions in the backward hemisphere 
w.r.t. the direction of the primary beam. A plausible model is one 
where the primary beam is directed downwards into the Sun, as in the 
thick-target model discussed above. 

As can be seen from Figures 4 and 5, F(He a /H 3 ) for the backward 
hemisphere can become as large as observed for He 3 -rich events. However, 
because this calculated enrichment occurs mainly at the very low energies 
^10 MeV/nucleon) , the isotopes have to be accelerated after their 
production. The need for post-production acceleration also follows 
from the fact that the He 1 * nuclei at the observed energies have much 
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shorter stopping ranges than the path length required for their 
production (Anglin, Dietrich and Simpson 1973b). Alternatively, 
in the thick-target model, the He 3 nuclei emitted in the backward 
direction probably have to go through a path length equal to the 
stopping range of the primary downward moving protons. Because these 
He 3 nuclei are produced by protons in the range 50 to 200 MeV, this 
path length is about 1 to 10 gcm -s . Since this range of values is 
much larger than the stopping length of the low energy He 3 nuclei, 
post production acceleration is again required. 

There is no data on the laboratory angular distribution of H . 
According to our kinematical assumption (Section II) , the tritium 
should be emitted in the forward direction; hence, as deuterium, it 
should not be observable in He 3 -rich events. However, even if some 
H 3 is emitted in the backward direction, it will be destroyed in nuclear 
reactions. 

Figure 12 shows the destruction cross sections of H 3 , H 3 and He 3 . 
Cross sections for the reaction H 3 (p,n)He 3 are from Hanson, Taschek and 
Williams (1949) and from Goldberg et al. (1961); cross sections for the 
reactions p+H 2 ->2p+n and p+He 3 -»d+2p are from Meyer (1972); and cross 
sections for the reaction He 3 (a?,p)Li s are obtained by using the detailed 
balance principle and data on the inverse reaction Li 6 (p,Q')He 3 (Yeronymo, 
Mani, and Sadeghi 1963; Marion, Weber and Mozer 1956; Hub, Clement and 
Wildermuth 1972) . As can be seen, H 3 has the largest destruction cross 
section at low energies; moreover, H 3 is transformed into He 3 . The 
destruction path length of H 3 is about 5 gem -3 , comparable to the range 
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traversed by the isotopes during post-production acceleration. We 
expect, therefore, that a large fraction of the backward moving 
tritium is transformed into He 3 . At low energies. He 3 can be 
destroyed only by alpha particles, but the weighted cross section is 
small . 

Finally, let us estimate the flux in the 2.2 MeV line that is 
expected from a He 3 -rich flare. Because both He 3 nuclei and neutrons 
are produced primarily in pet reactions (for primary spectra of the type 
observed from solar flares), we expect the production of about 1 neutron 
per He 3 nucleus. The efficiency of photon production is about 0.1 
(Wang and Ramaty 1974). About 30% of the He 3 is emitted in the backward 
direction (Section II). Furthermore, the escape probability of charged 
isotopes should be less than unity even if there is post-production 
acceleration. Therefore, we define a quantity f, the number of 2.2 
MeV photons per escaping He 3 nucleus, and allow it to be a free parameter 
of order unity. The time-integrated flux of 2.2 MeV photon; at earth 
is therefore 

/(2.2 MeV) = (4nR s ) -1 N (He 3 ) f (fih g g pg . ) , (20 ) 

where N(He 3 ) is the total number of He 3 nuclei released by the flare 
into the interplanetary medium. We estimate N(He 3 ) as follows: For both 

He 3 -rich events, October 14, 1969 and July 30, 1970, the alpha particle 
fluxes at about 10 MeV/nucleon were ~ 2 x 10 ~ 3 particles cm^s^sr’ 1 
(MeV/nucleon) 1 (M. Van Hollebeke, private communication)* Using a He 3 / 

He 4 ratio of about 0.3, we get that the He 3 flux at this energy is 
~ 6 x 10 4 . Then using the He 3 spectrum as given by Anglin (1974), we 
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get a total He 3 Intensity of ~ 10 “ 2 He 3 cm -3 s -1 sr 1 . This He 3 
flux is spread over some volume in the interplanetary medium. 

Ramaty and Lingenfelter (1973) have estimated this volume for the 
August 1972 events as ~ 10 39 cm 3 . Because the He 3 - rich events are 
very much smaller, we take the volume, somewhat arbitrarily, to be 
~ 5 x 10 37 cm 3 . Therefore, the predicted 2.2 MeV gamma-ray flux is 

^(2.2 MeV) “3 f N(He 3 )/10 3a ( £ ~* Ei ) • ( 21 ) 

This flux should be compared with the observed 2.2 MeV flux from the 
August 4, 1972 flare which was about 200 photons cm -3 . We hope that 
with improved instrumentation! nuclear gamma rays will be observed from 
He 3 -rich flares. 
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Table 1 

Isotope Production in pp,po? and op Reactions 


Deuterium Production 

E ^ (MeV/nucleon) 

E s (min) (MeV/nucleon) 

i. 

p + o* -> H 3 + 

23 


2. 

Qf + p -» H 2 + He 3 

23 

9.1 

3. 

p + a -* H 3 + p+n 

31.7 


4. 

Qf + p -* H 3 + p+n 

31.7 

12.5 

5. 

p + a ■* 2H 3 + p 

29.9 


6. 

a + p -* 2H 3 + p 

29.9 

11.8 

7. 

p + p -» H 3 + tt + 

286.5 

63.9 

Tritium Production 



8. 

p + a -» H 3 + 2p 

24.8 


9. 

a + p -» H 3 + 2p 

24.8 

13.1 

Helium-3 Production 



10. 

p + a -» He 3 + H 3 

23 


11. 

a + p -+ He 3 + H 3 

23 

12.2 

12. 

p + oi -+_He 3 + n + p 

24.8 


13. 

ot + p -* He 3 + n + p 

24.8 

13.1 
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Figure Captions 

1. Deuteron production in the thick-target model for exponential 

proton and alpha particle spectra with P^ = 150 MV and normal- 
ized to 1 proton of energy greater than 30 MeV. The numerals 
correspond to the reactions of Table 1, and backward indicates 
production into the entire backward hemisphere. 

2* Triton production in the thick-target model for exponential 

proton and alpha particle spectra with P Q - 150 MV and normal- 
ized to 1 proton of energy greater than 30 MeV. The numerals 
correspond to the reactions of Table 1. 

3. Helium-3 production in the thick-target model for exponential 

proton and alpha particle spectra with P^ - 150 MV and normal- 
ized to 1 proton of energy greater than 30 MeV, The numerals 
correspond to the reactions of Table 1* and backward indicates 
production into the entire backward hemisphere. 

4. The He 3 -to-H 3 ratio at the same kinetic energy per nucleon in the 

backward hemisphere calculated in the thick-target model with 
exponential spectra with various values of P Q . 

5. The He 3 -to-H 2 ratio at the same kinetic energy per nucleon in the 

backward hemisphere calculated in the thick-target model with 
power law spectra with various spectral indexes. 

6. Total deuterium production in the thin-target model for various 

primary spectra normalized to 1 proton of energy greater than 
30 MeV, Xj is the amount of matter traversed by relativistic 
particles . 
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7. Total tritium production in the thin-target model for various 

primary spectra normalized to 1 proton of energy greater 
than 30 MsV. is the amount of matter traversed by 
relavistic particles* 

8. Total helium-3 production in the thin- target model for various 

primary spectra normalized to 1 proton of energy greater than 
30 MeV. X l is the amount of matter traversed by relativistic 
particles • 

9. The He 3 -to-He 4 ratio at the same kinetic energy per nucleon in 

the thin-target model for various spectral indexes s- 

10. The H 3 -to-H x ratio at the same kinetic energy per nucleon in 

the thin-target model for various spectral indexes s. 

11. The H 3 -to-H 1 ratio at the same kinetic energy per nucleon in 

the thin-target model for various spectral indexes s. 

12. H 3 , H 3 and He 3 destruction cross sections. 
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